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Abstract 

A search for pair-production of long-lived, heavy, singly-charged particles has been 
performed with data collected by the ALEPH detector at a centre-of-mass energy of 
172 GeV. Data at ^/s = 161, 136, and 130 GeV are also included to improve the sensitivity 
to lower masses. No candidate is found in the data. A model-independent 95% conhdence 
level upper limit on the production cross section at 172 GeV of 0.2-0.4 pb is derived for 
masses between 45 and 86 GeV/C. This cross section limit implies, assuming the MSSM, 
a lower limit of 67 (69) GeV/C on the mass of right- (left-) handed long-lived scalar tans 
or scalar muons and of 86 GeV/C on the mass of long-lived charginos. 
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1 Introduction 


Most searches for supersymmetric (SUSY) particles assume that the lightest 
supersymmetric particle (LSP) is neutral and weakly interacting and that all charged 
SUSY particles ultimately decay, with a very small lifetime, into the LSP plus “standard” 
particles. In the hypothesis of R-parity conservation, this suggests missing energy as a 
possible signature of SUSY. 

Nevertheless, the possibility of supersymmetric long-lived charged particles which are 
not strongly interacting is favoured by some interesting classes of models. In Ref. [^, for 
example, a scenario is proposed in which the slepton is the next-to-lightest supersymmetric 
particle (NLSP). The slepton can decay into a standard lepton plus a Goldstino (£ —> IG) 
with a lifetime dependent on the SUSY-breaking energy scale \fF : 


cr ~ (130 yum) 


'100 GeV/c^ 


m-n 


5 


\ 

100 TeV } 
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Hence at LEP 2 centre-of-mass energies, if a/F is larger than a few thousand TeV, a 70 
GeV/c^ slepton has a decay length greater than 10 meters, which corresponds to the 
typical size of the detectors of the LEP experiments. 

Even in the minimal supersymmetric model (MSSM) a mass degeneracy of a few 
hundred MeV/c^ or less between the lightest chargino and the lightest neutralino could 
also induce a chargino lifetime long enough for the particle to decay outside the detector. 
This degeneracy may arise in some regions of the MSSM parameter space when relaxing 
the gauge unihcation condition (M = ^Mia.v?9w where M and M are the bino and 
wino masses). 

In R-parity violating SUSY models the possibility of long-lived charged weakly 
interacting particles also exists. If, for example, the slepton is the LSP a single dominant 
baryon-number violating coupling via the operator UDD would give a long lifetime to 
the slepton. 

At LEP 1, searches for heavy stable charged particles were performed by DELPHI |^, 
OPAL and ALEPH [Q, which sets the most stringent LEP 1 upper limit on the 
production cross section, roughly 1.5 pb at 95% GL for masses between 34 and 44 GeV/c^. 
DELPHI has analysed data collected at energies up to =172 GeV, setting a limit 
of ~ 0.3-0.5 pb for masses between 45 and 85 GeV/c^. 

This letter presents the results of a search for pair-production of long-lived, singly- 
charged, not strongly interacting particles with mass greater than 45 GeV/c^. The 
analysis is performed using data collected during the second 1996 LEP running period, at 
a/s = 172 GeV. To increase the sensitivity for low masses, data collected at 161 GeV (the 
first 1996 running period) and 130, 136 GeV (the 1995 high-energy run) are also analysed. 

The selection is based on kinematic criteria related to the pair-production hypothesis, 
and on the specihc energy-loss measurement {dE/dx), a powerful tool for exploring the 
high-mass region. When the mass of the heavy particles approaches the kinematic limit, 
the energy loss becomes high enough to saturate the main tracking detector’s electronics. 
To recover this interesting mass region, a selection based on the search for saturated 
signals in the tracking detector has been developed. 
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Limits on the cross section are translated into lower limits on the masses of long-lived 
slepton and charginos based on the prodnction cross sections predicted by the MSSM. 

2 The ALEPH detector 


A detailed description of the ALEPH detector can be fonnd in Ref. 0, and an acconnt of 
its performance as well as a description of the standard analysis algorithms in Ref. . 
Only a brief overview is given here. 

Charged particles are detected in the central part of the detector consisting of a silicon 
vertex detector, a drift chamber (ITC) and a time projection chamber (TPC), all immersed 
in a 1.5 T axial magnetic field provided by a snpercondncting solenoidal coil. A l/px 
resolution of 6 x 10“'^(GeV/c)“^ is measured. 

The TPC sense wires provide up to 338 measurements of the specihc ionization, dE/dx, 
for each track. To ensure a reliable dE/dx measurement, tracks in this analysis are 
required to have at least 50 associated wire hits. 

Between the TPC and the coil, an electromagnetic calorimeter (ECAL) is used to 
identify electrons and photons and to measure their energy, complemented by luminosity 
calorimeters (LCAL and SICAL) in the small polar angle region. The iron return yoke 
is instrumented to provide a measurement of the hadronic energy (HCAL) and, together 
with external chambers, muon identihcation. 

The two main ALEPH triggers relevant for this analysis are based on the coincidence 
between a track candidate in the ITC and an energy deposit in the ECAL or HCAL 
modules to which the track is pointing. 


3 Data and Monte Carlo samples 

Data collected at fys = 1720, 161, 136, and 130 GeV are analysed, corresponding to 
integrated luminosities of 10.6, 11.1, 2.9, and 2.9 pb“^, respectively. 

All major background reactions are generated at 161 and 172 GeV using the full 
detector simulation. These include the annihilation processes e’''e“ —>■ ffyy) and the 
various processes leading to four-fermion hnal states (e’''e“ ^ e’''e“ —>■ Wei/, 

e'''e“ —> Ze’''e“ and e’''e“ —>■ Zy*). These samples correspond to more than 100 times the 
integrated luminosity of the data. Two-photon processes (77 —^ and 77 —^ qq) are 

also simulated with an integrated luminosity about three times that of the data. 

For the signal, pair production of stable singly-charged particles is simulated for masses 
between 45 and 86 GeV/fy. These particles are treated as heavy muons by the detector 
simulation program. 


4 Low- and intermediate-mass selections 

A preselection is applied to reject topologies which are clearly incompatible with pair- 
production of back-to-back, massive, long-lived particles. Exactly two tracks are required, 

^The 172 GeV sample includes 1.1 pb“^ collected at 170 GeV. 
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both of which must be well reconstructed in the TPC. Both tracks must lie away from the 
beam axis (| cos6^1 < 0.90), have at least one hit in the ITC, and satisfy the conditions 
Moil + M02I < 0.3 cm and Moi| + M02I < 5 cm, where do (zq) is the distance of closest 
approach to the beam axis in the transverse plane (to the interaction point along the beam 
direction). The two tracks must have transverse momenta with respect to the beam axis, 
Pt, greater than O.ly^, equal momenta within three times the estimated error derived 
from the track fit, and the angle between them (acollinearity) must be greater than 160°. 
Both tracks must fail the electron identification cut |^, and deposit less than 20 GeV 
(50 GeV) in the EGAL (HGAL). An event is rejected if it contains a photon with energy 
above 250 MeV or if energy deposits are detected in the luminosity calorimeters. 

The acollinearity cut rejects radiative Z returns and two-photon processes. The 
cuts on Pt and the equal momentum requirement suppress two-photon and 
backgrounds. The photon veto rejects radiative Z returns and Bhabha events. The 
cuts on electromagnetic and hadronic energy further improve Bhabha rejection, the latter 
being relevant when the electrons enter an EGAL insensitive region. 

After this preselection, a large iJi~ background still survives, since the large mass of 
the signal has not yet been exploited. Two further sets of selection criteria, based on the 
measured particle masses mi, 2 , defined as m* = y^A'^eam “ Ph therefore introduced. 

The so-called low-mass selection requires the two particles to have masses in the range 
0.52 < mi ^/Ebearn < 0.80 and an acollinearity greater than 174°. In this mass range, 
the momenta of the tracks are high enough to produce an energy loss similar to that of 
ordinary particles; therefore dE/dx information is not used. 

The intermediate-mass selection requires masses in the range 0.80 < mi, 2 /Ebeam < 0.98 
and a specihc ionization + R^ 2 ) > 10. Here R^ is an estimator calculated by 
comparing the measured dE/dx, I, to that expected for a muon (/^): R^ = [I — (Ip)) /cj, 
where ai is the expected resolution of the measurement. In order that a single track 
saturating the TPG electronics does not suffice for the event to be selected, this estimator 
is assigned a value of 5 when saturation prevents a direct dE/dx measurement. 

For events with both masses in the range 90-98% of the beam energy the equal 
momentum requirement is relaxed, since the expected background (mainly e^e“ —> r+r“) 
is lower in this mass region. 

The only backgrounds surviving these cuts originate from the di-lepton channels, with 
cross sections of 10.4 fb for and 1.2 fb for r+r“. This corresponds to an expected 

background of 0.3 events in the full data sample. 

5 High-mass selection 

The extension of the search to masses approaching the kinematic limit (m/Ebeam — 1) 
requires a somewhat different approach than the low- and intermediate-mass analyses. 
Particles with masses within a few percent of the kinematic limit ionise so heavily that 
the TPG digitising electronics saturate and neither high-resolution spatial coordinates 
nor reliable dE/dx measurements are available. Saturation occurs at approximately 25 
times minimum-ionization. Saturated hits are not used by the standard ALEPH tracking 
since they have a resolution much worse than typical, unsaturated data. Nevertheless, the 
presence of many saturated channels is itself a distinctive feature of a high-mass signal 
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and in a clean two-track topology the characteristic pattern of a charged particle helix 
can be easily recognised from the spatial distribntion of saturated hits. 

The ITC amplifiers do not saturate, but the limited resolution in the z coordinate 
and the short lever-arm afforded by the ITC do not allow a kinematic selection based on 
tracks reconstructed with this detector alone. The high-mass selection requires exactly 
two tracks reconstructed in the ITC, which may or may not have (unsaturated) TPC hits 
associated to them. Events with energy deposited in the luminosity detectors are rejected. 
Energy deposits in the ECAL and HCAL are used to provide a three-dimensional point on 
each track candidate and verify a roughly back-to-back topology. The two most energetic 
calorimeter objects are selected; they must have an acollinearity of at least 160°, opposite 
polar angles within 8.5°, and each contain at least 50% of the calorimetric energy deposited 
on their side of the detector (dehned by a plane containing the interaction point and 
perpendicular to the beam direction). Since low-/?, massive particles should not shower 
in the calorimeters, the energy of these objects must be less than 5 GeV. 

The positions of the two selected calorimeter objects, used as estimates of where the 
particles may have exited the tracking volume, and the position of the interaction point, 
used as an estimate of the origin of these particles, provide enough information for a single- 
helix fit of the transverse momentum and polar angle of the two track candidates. The 
reconstructed helix must cross at least four TPC pad rows on each side of the detector. 
The particle masses mi 2 must be between 90 and 99.5% of the beam energy. 

The estimated particle trajectories from the single-helix fit do not rely on any 
information from the tracking detectors, so the pattern of all hits (saturated and 
unsaturated) in the ITC and TPC can be checked for consistency. An approximatively 
5 standard deviations road is dehned around the single-helix trajectory; this road has a 
width of 1 cm (in rcfi) for the ITC, and 5 cm (in both ref) and z) for the TPC. At least 
60% of the ITC layers and 40% of the TPC pad rows crossed by the road must contain a 
hit; this conhrms the presence of charged tracks where predicted by the ht to calorimeter 
objects. The fairly low occupancy requirement in the TPC allows a track which projects 
onto a crack between outer sectors to be retained if every pad row in the inner sectors 
contains a hit. 

Finally, the high mass of both particles is conhrmed by the dE/dx measurements 
(if available) or by saturation of the TPC wires or pads (if not) by requiring either 
-|- R ^2 > 10 or /sat > 0.33; /sat is the fraction of TPC pad rows crossed by the htted 
helix which contain a saturated coordinate within the 5 cm road. 

No background from the simulation survives the high-mass selection. 


6 Efficiency 

The hnal selection is a combination of the low- and intermediate-mass criteria (including 
their common preselection), and the complementary high-mass analysis. Together, the 
three selections yield an efficiency rising gently from 50% at m/Ebeam = 0.50 to 70% 
at m/Ebearn = 0.90; this is shown in Fig. |I| where the selection efficiency (including the 
trigger) for spin-1/2 stable charged particles is plotted as a function of m/Ebeam- 

For m/Ebearn > 0.993 the particles stop in the calorimeters because of the high energy 
loss. In this case, if the particles are not stable and decay in the calorimeters, they 
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might produce additional energy deposits which will affect the selection efficiency. In 
the Goldstino scenario (long-lived sleptons) the energy release can be important and the 
efficiency is conservatively assumed to be zero. However in the high mass-degeneracy 
scenario (long-lived charginos) almost all the energy is taken by the neutralino and only a 
few hundred MeV are deposited in the calorimeters. In this case the efficiency is assnmed 
to be eqnal to that of stable particles. 

The efficiency depends on the spin (i.e., the production angular distribution) and 
mass of the charged particle. In the Monte Carlo signal generation, massive particles are 
prodnced with a flat angnlar distribntion. The spin-1/2 and spin-0 efficiencies are derived 
from the generated signal by rescaling the angular acceptance using the formulae 


da (spin = 0) 

dC 

da(spin = 1/2) 

dC 


oc 

oc 



s 


l + cos^9 + {l-p^)sm^e 


( 1 ) 


The expression for the differential spin-1/2 cross section is valid for s channel production 
under the assumption of purely vector coupling to spin-1 bosons. 

To acconnt for the dependence on ^/s, all cuts are performed in terms of variables 
normalised to the beam energy and the efficiency is parametrized as a fnnction of m/i^beam- 
The trigger efficiency has been estimated with a Monte Carlo program simulating the 
two ALEPH triggers briefly described in Section 2. The trigger efficiency exceeds 99% 
for m/Ebearn < 0.97, crosses 99% at m/Ebeam — 0.98, and finally decreases to 80% for 
"i/^beam > 0.99. 


7 Systematic checks 


For all variables used in the selections, good agreement exists between data and Monte 
Carlo simnlation. The most important variables used to isolate the signal events are the 
measnred momenta and dE/dx of the two particles. The nncertainty on the beam energy 
has a negligible impact on the results. The tracking performance was tested with data at 
the Z peak with a variety of techniques (see for example |^), showing agreement of the 
absolute momentum calibration between data and Monte Carlo better than ~ 100 MeV/c 
for 45 CeV/c tracks. This level of momentnm nncertainty has a negligible impact on the 
analysis. 

Stndies of dE/dx performance at the Z peak have been described in many ALEPH 
papers (see for example 0). Additional studies were performed with the 1996 data, 
using minimnm-ionising pions, identified muons, and electrons identified using the ECAL 
information. Agreement of the absolute dE/dx calibration is obtained between data and 
Monte Carlo at the level of ~ 0.3 times the expected resolution, which is ~ 4.5% for high 
energy isolated electrons. The 1996 data exhibit a dE/dx resolution 10-20% better than 
the Monte Carlo which goes in the direction of increasing the selection efficiency. These 
effects lead to a very small change (< 1%) of the selection efficiency since the dE/dx 
information is only used in a /? region when the expected ionization is much larger than 
the cut applied in the analysis. 
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For the high-mass selection, the identihcation of pair-produced charged tracks using 
calorimeter objects has been tested using di-muon events collected at the Z peak. 
Although muons produced in Z decays do not ionise heavily enough to saturate the TPC 
electronics, the high-mass signal can be simulated with data by applying to di-muon hits 
the same algorithm used for saturated coordinates. The resulting TPC hit positions agree 
with those of the standard ALEPH reconstruction to within about half the width of a 
cluster on the pad row in ref) (~ 1 cm) or half the pulse-length in z (-C 1 cm). 

The transverse spatial resolution of the ECAL entry point used in the fit is roughly 
1.5 cm, consistent with the granularity of the ECAL towers. The longitudinal uncertainty 
on the entry point is expected to be of the order of one radiation length in lead, or less 
than 2 cm, and in any case has only a second-order effect on the ht. 

A conservative 5% systematic error on the selection efficiency has been applied. The 
impact of this systematic error on the cross section upper limit has been calculated 
according to the method described in Ref. [^, based on the convolution of the Poisson 
probability estimator with a Gaussian of sigma equal to the error on the signal efficiency. 

8 Results 

No event in the 172 GeV data survives any of the three selections. This sets the 95% 
confidence level upper limits on the production cross section which are plotted as dashed 
curves in Fig. ^ and ^ for spin-1/2 and spin-0 particles, respectively. No candidate is 
found in the data samples collected at ^/s = 161, 136 and 130 GeV, allowing the cross 
section limit to be improved. 

The corresponding integrated luminosities must be rescaled to account for the 
dependence on ^/s of the spin-0 and spin-1/2 production cross sections shown in Eq. ^ 
The combined upper limit at 95% GL is 

3 

{^, 172^172 + + Afe^ise + Afo^iso) 

where denotes the rescaled luminosities and e the selection efficiencies. The improved 
limits obtained by including the data at lower energies are shown as thick curves in Fig. ^ 
and ^ Gross sections of 0.2-0.4 pb at ^/s = 172 GeV are excluded at 95% GL for masses 
between 45 and 86 GeV/c^, almost independently of the mass and spin of the heavy 
particle. 

As already mentioned, the possibility of a long-lived charged scalar lepton depends, 
in the model of Ref. |^, on the value of the SUSY breaking scale parameter a/F; the 
requirement that the particle decay outside the detector therefore constrains this scale. 
In Fig. ^ the cross section limit is plotted as a function of m and a/F. The lifetime- 
dependence of the efficiency is approximated with a factor exp (—2f'det//^qcr), where 
^det = 8.5 m is the maximum length travelled by a particle inside the detector. The 
limit on the cross section worsens for scales below 500-1000 TeV. 

These cross section upper limits can be translated into lower limits on the masses of 
charginos and sleptons with long lifetimes. The cross sections expected in the MSSM 
for charginos and for right- and left-handed smuons or staus are also shown in Fig. ^ and ^ 
respectively. The selectron production cross section also depends on the neutralino masses 
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and couplings due to neutralino exchange in the t channel, therefore no general mass limits 
can be given. For charginos the cross section range displayed in Fig. |]was determined as 
follows: a value of tan [3 = \/2 was chosen; a given chargino mass then dehnes a relation 
between M and /i (where fi is the supersymmetric mass term which mixes the two Higgs 
superfield); for each pair of M and /i values, the gauge unification condition was relaxed 
and M varied up to 5 M; if the condition — m^o < 200 MeV/c^ could be satisfied 
in that way, the chargino production cross section was calculated, assuming a sneutrino 
mass of 250 GeV/c^. 

The lower 95% confidence level mass limits for long-lived smuons or staus are 67 and 
69 GeV/c^ for right- and left-handed particles, respectively. For a long-lived chargino, due 
to the large production cross section, the kinematic limit of 86 GeV/c^ is almost attained. 

9 Conclusions 

ALEPH data collected at 172, 161, 136 and 130 GeV yield a 95% conhdence level upper 
limit of 0.2-0.4 pb on the pair-production cross section at 172 GeV of long-lived, singly- 
charged particles with masses between 45 and 86 GeV/c^. This cross section limit implies, 
in the MSSM, lower limits of 67 GeV/c^ on the mass of right-handed staus and smuons 
and of 69 GeV/c^ on the mass of left-handed staus and smuons. For long-lived charginos 
a lower mass limit of ~ 86 GeV/c^ has been obtained. 
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Figure 1: Selection efficiency as a function of for spin-1/2 stable charged 

particles. The efficiency curve for spin-0 particles has a similar shape but is slightly 
higher, owing to the larger fraction of cross section within the geometrical acceptance 
of the selection. 
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Figure 2: 95% CL upper limit on the production cross section for spin-1/2 particles at 
a/s = 172 GeV, as a function of mass. The dashed curve shows the limit obtained from 
the 172 GeV data set alone, while the thick solid curve shows the limit obtained by adding 
the 161, 136 and 130 GeV data to the analysis. The band on the right shows the chargino 
cross section calculated as described in Section 8. 
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Figure 3: 95% CL upper limit on the production cross section for spin-0 particles at 
a/s = 172 GeV, as a function of mass. The dashed curve shows the limit obtained from 
the 172 GeV data set alone, while the thick solid curve shows the limit obtained by 
adding the 161, 136 and 130 GeV data to the analysis. The left-handed and right-handed 
scalar-muon (or scalar-tau) production cross sections are superimposed. 
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Figure 4: 95% CL upper limits (in pb) on the production cross section for spin-0 particles 
at a/s = 172 GeV as a function of mass and SUSY-breaking scale, a/F. 
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